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The classical Heisenberg antiferromagnet with uniaxial exchange anisotropy, the XXZ model, and 
competing planar single-ion anisotropy in a magnetic field on a simple cubic lattice is studied with 
the help of extensive Monte Carlo simulations. The biconical (supersolid) phase, bordering the 
antiferromagnetic and spin-flop phases, is found to become thermally unstable well below the onset 
of the disordered, paramagnetic phase, leading to interesting multicritical points. 
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Uniaxially anisotropic Heisenberg antiferromagnets in 
a magnetic field have been studied quite extensively in 
the past, both experimentally and theoretically [l[ 
Usually, they display, at low temperatures and fields, the 
antiferromagnetic (AF) phase and, when increasing the 
field, the spin-flop (SF) phase. A prototypical model de- 
scribing these phases as well as multicritical points, is the 
Heisenberg model with a uniaxial exchange anisotropy, 
the XXZ model 

n xxz = jJ2[&(s?s? + SfS*) + SfS?] - H sf 

(1) 

where J > is the antiferromagnetic exchange coupling 
between classical spins, (^?(j)>^f(j)>5f(j)); of length one 
at neighboring sites, i and j, on a simple cubic lattice, A 
is the uniaxial exchange anisotropy, 1 > A > 0, and H 
is the applied magnetic field along the easy axis, the z— 
axis.- Recently, we presented numerical evidence for the 
existence of a bicritical point in the Heisenberg universal- 
ity class Q, at which the AF, SF and paramagnetic (P) 
phases meet. Especially, there the value of the critical 
isotropic Binder cumulant agrees closely with the one of 
the perfect Heisenberg model. The scenario is in accor- 
dance with early renormalization group (RG) analyzes in 
one-loop-order [1,0], as well as a subsequent RG calcula- 
tion in two-loop-order [(| showing that such a bicritical 
point may exist. Note that a conflicting characteriza- 
tion of the multicritical AF-SF-P point as a 'biconical 
fixed point' @, [1[ has been suggested, mainly based on 
high-order RG arguments, augmented by Monte Carlo 
simulations. 

In this contribution, we shall consider the XXZ model, 
Eq. (1), with an additional planar single-ion anisotropy, 
competing with the uniaxial exchange anisotropy. This 
term reads 
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FIG. 1. Spin orientations on neighboring sites showing an- 
tiferromagnetic (a), biconical (b), and spin-flop (c) ground 
state configurations. 

magnets. One may also mention a study for the clas- 
sical full model on square lattices [l4j. Again, the BC 
phase arises from the competition between the uniaxial 
exchange and planar single-ion anistropies. 

In the BC phase, see Fig.l, spins on the bipartite sub- 
lattices form cones with different opening angles so that 
one gets coexisting orders of the staggered transverse, 
xy, spin components, like in the SF phase, as well as 
longitudinal, z, components, like in the AF phase. 

The resulting ground state phase diagram for the clas- 
sical model on a simple cubic lattice may be determined 
analytically [IH, [l4j . Possible ground states are depicted 
in Fig. 1, showing AF, SF, and BC spin configurations. 
A concrete example of the ground state phase diagram 
in the (H/J, D/J) plane, setting A = 0.8, is displayed in 
Fig. 2. The BC structure is seen to intervene between 
the AF and SF structures for non-zero planar single-ion 
anisotropy D/J < 0.6. The opening angles ®a and Qb 
are uniquely related, depending on the field, continuously 
approaching the limiting values of the AF and SF orien- 
tations. The BC structure reaches its maximal extent 
close to D/J = 0.4. The lower critical field between the 
AF and BC ground states, H ab , is given by 



U si = Dj2(Sif (2) 

i 

The full model, T-i xxz + ~risi , has been investigated, in 
the last few years, carefully, especially, for S = 1 quan- 
tum antiferro mag netic spin chains, displaying a super- 
solid phase [ol— Fill] . being the analog of the mixed or bi- 
conical (BC) phase 0, [12, HH in standard notation for 



H ab = V(6J-2£) 2 -(6JA) 2 (3) 

and the upper critical field, separating the BC and SF 
ground states, H bs , is given by 

H bs = [ 36 J 2 - (6 J A + 2D) 2 ] /H ab . (4) 
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FIG. 2. Ground states in the (H/J,D/J) plane of the full 
model on a cubic lattice with uniaxial exchange anisotropy 
A = 0.8 and planar single-ion anisotropy, D/J > 0. 

The transition from the SF to the ferromagnetic (F) 
ground state occurs at the field H s f with 

H sf = 6 J(l + A) + 2D. (5) 

At non- vanishing temperatures, the BC phase may be 
expected to lead to a tetracritical point [l5|, at which 
the AF, SF, BC, and P phases meet. One of the aims of 
the present study has been to clarify this expectation by 
using extensively Monte Carlo (MC) techniques (l6| . 

Followin g pr evious MC studies on the XXZ model and 
variants [1,13, Qjl [l8| , we fixed the exchange anisotropy 
at A = 0.8. The single-ion anisotropy has been chosen to 
be, mostly, D/J = 0.2. At the end of this paper, results 
for D/J = 0.4 will be briefly discussed as well. To take 
into account systematically finite-size effects [l9| , we sim- 
ulated cubic lattices with L 3 spins, with L ranging from 
4 to 32, employing full periodic boundary conditions. To 
obtain data of the desired accuracy, runs with, at least, 
5x 10 7 Monte Carlo steps per spin were performed, apply- 
ing the standard Metropolis MC algorithm. Error bars 
result from averaging over a few of such runs. When er- 
ror bars are not depicted in the figures, then they are 
smaller than symbol sizes. 

To determine the phase diagram, we record, especially, 
quantities related to the Ising, XY, Heisenberg, and bi- 
conical order parameters. In particular, we monitor the 
absolute values of the staggered longitudinal, m z t , trans- 
verse, m s t , and isotropic, m x J! z , magnetizations as well as 
of the product of the staggered longitudinal and trans- 
verse magnetization components. Furthermore, related 
staggered susceptibilities, namely, Xsf X°st ' 

and XsV i are 

computed. In addition, we calculate corresponding longi- 
tudinal, transverse, and isotropic Binder cumulants [20j, 
U z , U xy , and U xyz . 

The interesting part of the resulting phase diagram 
is depicted in Fig. 3: Two distinct multicritical points 
occur, the AF-SP-P and the AF-BC-SF points. 







1 1 1 
















* 




3.5 










! P 






SF 




at 




3,4 
























=e 












3,3 






Jg'' 






3,2 






AF 






3.1 c 




1,1 








0,4 0,5 


0,6 0,7 


0,8 0,9 


i 


1,1 1,2 



k n T/J 



FIG. 3. Phase diagram in the temperature-field plane of the 
three-dimensional full model with A = 0.8 and planar single- 
ion anisotropy D/J — 0.2. The broken line indicates transi- 
tions of first order. 

As before Q, to obtain the AF-P boundary, with the 
transition being in the Ising universality class, the size 
dependence of the height and position of the maximum 
in the longitudinal staggered susceptibility, Xsti allows 
for reliable estimates. U xy turns out to be very useful in 
estimating the SF-P transition line. There one observes 
rather small finite-size correction terms to the critical 
Binder cumulant for cubic lattices in the XY universal- 
ity class [U, (U xy )* = 0.586.., with the critical Binder 
cumulant being the cumulant at the transition in the 
thermodynamic limit, L — » oo. It turned out to be less 
straightforward to identify the transition between the AF 
and SF phases, with, possibly, an intervening BC phase. 

Indeed, the extent of the BC phase follows from the 
behavior of the staggered longitudinal and tranverse sus- 
ceptibilities, Xst anc l Xst i w ith the AF-SF transition be- 
ing signalled by the peak in x*t> anc l the BC-SF tran- 
sition being indicated by the location of the maxima in 
the longitudinal susceptibility. Obviously, the BC phase 
shrinks with increasing temperature, becoming, even- 
tually, very narrow. Accordingly, further analyzes are 
needed to check the possible vanishing of the BC phase, 
which would then imply the existence of an AF-BC-SF 
multicritical point. 

A convenient way for identifying the AF-BC-SF point 
is demonstrated in Fig. 4. There, the size dependence of 
the height of the maxima in the staggered longitudinal 
and transverse susceptibilities is shown, both in the case 
of an intervening BC phase, at the lower temperature, 
ksT/ J = 0.6, and in case of a direct AF-SF transition of 
first order, at the higher temperature, ksT / J = 0.9. 

In fact, at the higher temperature, the maxima of both 
susceptibilties grow to a good approximation like L 3 , for 
L exceeding about 8, see Fig. 4. This behavior is in 
accordance with a transition of first order. In marked 
contrast, at the lower temperature, ksT / J = 0.6, both 
maxima occur at distinct fields, reflecting the separate 
AF-BC and BC-SF transitions. Moreover, the heights of 
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FIG. 4. Finite size behavior of the maximal staggered longi- 
tudinal, z (open symbols), and transverse, xy (full symbols), 
susceptibilities for linear dimension, L, ranging from 4 to 24, 
at ksT/J = 0.6 (circles) and 0.9 (squares), with A = 0.8 and 
D/J — 0.2. For comparison, the dependence oc L 3 , charac- 
teristic for a transition of first order, is shown as solid line. 




ii i i i i i i i i i i i i 

.,45 3,46 3,47 3,48 3,49 3.5 3,51 
H/J 



FIG. 5. Isotropic Binder cumulant U xyz , for A = 0.8 and 
D/J = 0.2 at k B T/J = 1.075, varying the field, for lattice 
sizes L = 4 (circles), 8 (squares), 12 (diamonds), 16 (triangles 
up), 20 (triangles left), 24 (triangles down), and 32 (triangles 
right). 

the peaks grow much weaker with system size: The effec- 
tive exponents, described by the local slopes between suc- 
cessive points in the doubly logarithmic plot, are clearly 
smaller than 3, as illustrated in Fig. 4. Actually, for both 
susceptibilities we find similar results for the effective ex- 
ponents, dropping down, for the largest sizes we studied, 
to values even well below 2, to about 1.6. 

Actually, one may expect the AF-BC transition to 
belong to the XY, and the BC-SF transition to the 
Ising universality class 15!]. A detailed description of 
these transitions is, however, beyond the scope of the 
present study. We also refrained from attempting to lo- 
cate the AF-BC-SF point very accurately, and to eluci- 
date its (critical) properties. Its position is roughly at 
k B T abs /J = 0.85 ±0.05. 

In any event, the AF-BC-SF point occurs well below 



the AF-SF-P multicritical point, k B T asp /J = 1.075 ± 
0.01, H asp /J = 3.490 ± 0.002, at which the continuous 
AF-P and SF-P phase transition lines meet the AF-SF 
transition line of first order. Obviously, the upper mul- 
ticritical point reminds of the multicritical point in the 
XXZ model on the cubic lattice. Accordingly, the planar 
single-ion anisotropy leads to a BC phase, but it does not 
lead to a tetracritical AF-BC-SF-P point, at least, in the 
case we considered so far, D/J = 0.2. 

To elucidate the nature of the AF-SF-P multicritical 
point, we studied in its vicinity, especially, the isotropic 
Binder cumulant U xyz (T, H, L). In case of a bicriti- 
cal point with Heisenberg symmetry, the critical Binder 
cumulant {JJ xyz )* may be expected [12, HH to acquire 
the characteristic value 0.620 (±0.003), using periodic 
boundary conditions for lattices of cubic shape. 

Similarly to the findings for the perfect XXZ model, 
D = 0, we may distinguish three different scenarios, when 
varying the temperature and monitoring the size depen- 
dence of U xyz as a function of the field H [3] . 

Below the AF-SF-P point, e. g. , at k B T/J = 1.05, 
U xyz increases with larger lattice size L, so that there 
are no intersection points of the isotropic cumulants near 
the direct AF-SF transition. 

At k B T/J = 1.075, see Fig. 5, all intersection points 
of the isotropic cumulants, for lattices of sizes L ranging 
from 8 to 32, occur closely to the characteristic Heisen- 
berg value, (U xyz )* « 0.620. This fact may be inter- 
preted as evidence for being in the immediate vicinity 
of a bicritical point of Heisenberg symmetry. Actually, 
moving to higher temperatures, another scenario shows 
up. For example, at k B T/J = 1.1, there are still inter- 
section points of the cumulants for different lattice sizes, 
but they shift systematically for larger lattices to lower 
and lower values below that of the critical cumulant in 
the Heisenberg case. For instance, the intersection point 
of the cumulants for the two largest sizes, L = 24 and 32, 
is even below 0.58. 

As for the perfect XXZ model, these observations on 
the isotropic Binder cumulant are consistent with the 
AF-SF-P point, at D/J = 0.2, being a bicritical point 
in the Heisenberg universality class. 

Note that the staggered isotropic susceptibility seems 
to be the less sensitive quantity to identify the nature of 
the AF-SF-P point, with the values of the critical expo- 
nent for bicritical and biconical multicritical points being 
very close to each other [8j. 

Finally, we mention briefly results of our (preliminary) 
simulations for an enhanced planar single-ion anisotropy, 
D/J = 0.4. Again, there is no evidence for a tetracrit- 
ical AF-BC-SF-P point. Instead, we find an AF-BC-SF 
point. It is followed, at higher temperatures, by a direct 
transiton of first order between the AF and SF phases, 
with the BC phase being squeezed out. Upon further 
increase of the temperature, an AF-SF-P multicritical 
point occurs. 

I should like to thank R. Folk and M. Hasenbusch for 
useful discussions. 
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